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ABSTRACT 


The  theory  of  automatic  detection  at  a  seismic 
array  by  means  of  an  F  statistic  is  extended  to  a  net¬ 
work  of  such  arrays.  The  arrays  may  have  equal  or 
different  expected  signal-to-noise  ratios.  Two  tech¬ 
niques  are  discussed:  (1)  the  composite  F  detector  in 
which  a  vote  is  taken  among  the  arrays;  (2)  a  multi¬ 
array  F  detector  in  which  the  original  data  from  the 
independent  arrays  are  combined  to  form  one  F  statistic. 
The  detectors  are  found  to  be  nearly  equal  in  detection 
capability,  with  the  multi-array  detector  superior  by 
1-2  dB  in  the  cases  examined.  For  example,  with  22 
independent  arrays  of  6  elements  each,  assuming  equal 
expected  signal-to-noise  values,  the  two  detectors  are 
5.4  and  4.3  dB  worse,  respectively,  than  an  F  detector 
would  be  operating  on  a  beam  of  6x22-132  channels  with 
perfect  signal  correlation. 
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INTRODUCTION 


The  performance  of  an  F-detector  on  a  31-element 
TFO  short-period  array  has  been  studied  by  Blandford 
(1970).  lie  found  that  the  detector  operated  in  good 
agreement  with  theory;  therefore  it  seems  reasonable  to 
study  the  expected  performance  of  an  F-detector  operat¬ 
ing  on  a  network.  This  would  be  of  application  to  at 
least  two  distinct  types  of  systems:  (1)  an  array  such 
as  NORSAR  where  signal  correlation  between  subarrays 
may  be  poor;  in  this  case  a  detector  might  be  operated 
independently  for  each  subarray,  and  the  results 
combined  in  some  way;  (2)  a  world-wide  network  of 
completely  independent  arrays  with  different  noise  and 
signal  levels. 

A  simple  F-detector  is  said  to  detect  if  FONj,^,*) 
^  Fq,  where  F(N1,N2,A)  is  the  computed  test 
statistic,  assumed  to  belong  to  a  population  distributed 
as  a  non-central  F  with  and  N2  degrees  of  free¬ 
dom  and  noncentrality  parameter  A,  and  Fq  is  the  thres¬ 
hold  determined  from  the  desired  false  alarm  rate.  Thus 
the  probability  of  detection  by  this  detector  is  the 
probability  that  F  Fq  in  the  presence  of  a  signal.  In 
order  to  make  meaningful  comparisons  between  different 
systems,  the  false  alarm  rate  should  be  held  constant, 
and  this  is  done  in  all  the  measurements  made  in  this 
report.  Blandford  (1970)  suggests  that  =  2BT  =  3  is 
suitable,  where  integrated  bandwidth  =  B  =  0.5  Hz,  and 
the  signal  window  =  T  =  3.0  seconds.  Theory  shows  that 
N2  =  2BT (N-l),  where  N  is  the  number  of  elements  in  the 
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array.  (Sonc  uotational  confusion  is  possible  with  the 
denominator  of  the  expression  for  signai-tc-noisc  ratio 
(S/N);  however,  this  is  usually  clear  from  the  context.) 

A  power  detector  nay  he  defined  as  the  ratio  of 
the  beam  power  in  a  short  signal  window  to  that  in  an 
arbitrarily  long  noise  window.  If  the  noise  is  station¬ 
ary  this  is  equivalent,  for  small  signals,  to  at.  F 
detector  with  X,  arbitrarily  large,  approaching  a 
distribution.  (liven  the  hypothesis  of  stationary  noise, 
this  would  he  superior  in  pure  detection  capability  to 
the  l'-detcctor,  where  the  noise  statistics  are  accumu¬ 
lated  only  in  the  signal  window,  and  where  is  there¬ 
fore  limited.  However,  Blandford  (1970)  showed  that  the 
I  detector  lias  features  which  make  it  more  reliable  in 
practice,  and  in  this  paper  we  show  that  by  use  of  long 
time  windows  the  quantitative  difference  can  be  made 
negligible  for  all  except  the  smallest  arrays  of  2-3 
elements . 

Signal-to-noisc  values  used  throughout  this  report 
are  calculated  as  the  square  root  of  power  ratios. 
Blandford  showed  that  to  convert  these  values  to  those 
given  by  the  ratio  of  pcak-to-pcak  signal  divided  by 
peak-to-pcak  noise  in  the  previous  60  seconds,  the 
square  root  of  the  power  ratio  should  be  multiplied  by 
0,5.  Signal-to-noisc  values  are  defined  at  both  the 
instrument  and  beam  level  by  the  symbols  (S/N)scjs  and 
(f/^boam  respectively.  The  relation  fS/M)bcam  = 

\l '“(S/N)  .  is  true  since  the  noise  is  uncorrelatcd , 

hv  hypothesis  for  the  F-dctector. 
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In  the  following  sections,  we  first  discuss  the 
question  of  the  optimum  signal  window  length  over  which 
to  calculate  the  P  statistic.  Blandford  used  3  seconds, 
hut  this  is  shown  not  to  he  optimum  for  small  arrays. 

he  then  develop  the  theory  for  a  “composite  I  - 
detector"  in  which  a  vote  is  taken  among  independently 
detecting  arrays  in  order  to  settle  on  a  final  detec¬ 
tion. 

This  is  followed  hy  a  section  in  which  the  "multi¬ 
array  P-dctcctor"  is  developed.  By  this  technique  the 
actual  traces  from  each  independent  array  are  combined 
into  a  single  statistic  which  for  uncorrclatcd  signals 
is,  in  theory,  superior  to  any  other  detector,  includ¬ 
ing  the  composite  detector. 

The  two  methods  arc  then  extended  to  the  ease  of 
unequal  expected  signal-to-noisc  ratios  at  the  diff¬ 
erent  subarrays.  This  is,  of  course,  particularly 
applicable  to  a  world-wide  network  of  stations  some  of 
which  would  be  closer  to  an.  epicenter  than  others. 
However,  many  workers  have  also  shown  that  for  small 
seismic  regions  of  the  earth  there  arc  significant 
repeatable  amplitude  anomalies  recorded  across  I.ASA, 
so  the  two  methods  would  also  be  suitable  detectors 
for  such  a  case. 
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01'  I  I  MUM  SICNAt.  ft  I  NDOh  1  OR  IHiTI  CTlON 


f landlord  (1970)  used  a  3-second  signal  window  in 

liis  study  of  the  automatic  detector  at  Tl'O.  This  window 

w  is  suggested  by  the  studies  of  earlier  workers  and  by 

intuition,  which  suggests  that  the  window  should  be 

only  as  long  as  the  dominant  portion  of  the  signal.  In 

Ins  figure  S,  Itlandford  showed  that  for  small  arrays, 

\  7,  there  were  severe  departures  from  proportionality 

1  /  ’ 

to  S  for  the  s ignal-to-noisc  ratio  which  would  be 
detected  DO  per  cent  of  the  time. 

It  seemed  plausible  that  this  departure  occurred 
because  the  noise  estimate,  calculated  in  the  signal 
window  by  averaging  residual  power  over  all  N  channels, 
was  not  statistically  reliable  for  small  X.  The  sta¬ 
bility  would  he  increased  by  a  longer  time  window,  but 
the  s ignal-to-noisc  ratio  would  decrease  by  (3/T)1^. 

It  turns  out,  as  we  shall  see  below,  that  the  two 
effects  almost  perfectly  cancel  for  large  arrays,  and 
that  detection  is  substantially  improved  for  small 
arrays . 

y 

As  T  increases,  \  »  2BT (S/N) remains  constant, 
since  (S/X)“  -  1/T.  X,  ami  X,  increase  in  proportion 
to  T/3.  Calculation  for  successive  values  of  T  yields 
the  curves  in  figure  i.  In  the  neighborhood  of  the 
optimum  window  length  the  threshold  signal  is  insensi¬ 
tive  to  the  precise  value  of  the  window  length.  However, 
for  3  elements  there  is  almost  0.5  m^  improvement  in 
the  required  signal  between  the  3-sccond  and  30-sccond 
w i udows . 


-I 


The  optimum  window  length  us  u  function  of  the 
number  of  elements  is  shown  in  figure  2.  The  actual 
optimum  length  is  presumably  somewhat  longer  since  in 
reality  the  signal  is  not  exactly  zero  outside  of  a 
3-sctond  window. 

In  Figure  3  the  upper  curve  gives  the  element 
signal/noise  value  required  for  90  per  cent  prob¬ 
ability  of  detection  at  0.1  false-aiarm  per  day  for  a 
3  second  time  window  as  a  function  of  the  number  of 
elements  in  the  array.  The  lower  curve  gives  the 
signul/noine  value  for  the  optimum  window.  The  curves 
arc  shifted  up  by  a  factor  of  /I  over  those  given  h> 
((landlord  (1970)  due  to  an  error  in  calculation  in 
that  paper. 

1  / 

We  note  that  both  curves  arc  asymptotic  to  N 
the  theoretical  performance  curve  for  a  power  detector. 
This  can  he  simply  explained  by  noting  that  the  signal/ 
noise  ratio  required  on  the  beam  for  a  given  probabi  lity 
of  detection  decreases  as  a  function  of  N,  but 
approaches  a  limit  for  large  N,  for  which  further 
increases  of  ^  have  little  effect.  Thus  for  large  N , 
the  decrease  of  seismometer  (S/N)  required  results 
primarily  from  the  improvement  of  (S/N)  on  the 

beam. 
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AT  LkAST  K  SUBARRAY  DI-TLCT  ION 


Mien  the  signal  is  not  identical  across  a  large 
array,  as  may  he  the  case  at  VORSAR,  the  array  can  he 
broken  into  a  number  of  smaller  suharrays  and  the  I- 
detector  operated  separately  on  each  subarray.  If  one 
then  asks  for  at  least  K  suharrays  to  detect  before  one 
declares  a  detection  (K  is  a  number  yet  to  he  specified) 
then  the  detection  capability  is  not  much  worse  than 
one  would  obtain  in  the  ideal  case  of  an  identical 
signal  on  the  entire  array  acting  as  one  !•'  detector. 

It  turns  out  that  there  is  an  optimum  choice  of  K  which 
maximizes  the  signal  detection  probability  with  fixed 
false-alarm  rate.  This  value  is  usually  slightly  more 
than  half  the  number  of  suharrays,  and  it  depends  on 
the  numbers  of  degrees  of  freedom  (i.c.,  on  the  time 
and  frequency  windows  and  the  number  of  elements  per 
subarray)  and  on  the  number  of  subarrays,  also  to  a 
slight  extent  on  the  chosen  false-alarm  rate  (F  thres¬ 
hold)  and  signal  level.  The  maximum  is,  however, 
fairly  broad  and  an  error  of  1  or  2  in  the  choice  of 
k  would  not  seriously  degrade  the  performance  of  a 
large  system, 

he  assume  in  our  analysis  that  the  signals  arc 
perfectly  correlated  across  each  individual  subarray 
(and  the  noise  perfectly  uncorrclatcd) ,  In  addition, 
we  will  begin  by  assuming  the  signals  on  different 
suharrays  to  he  of  the  same  sire  (same  signal/noise 
ratio),  although  later  we  will  consider  the  effects  of 
different  sice  signals.  Representing  the  probability  of 
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a  detection  on  any  subarray  by  p  (the  same  for  all  sub- 
arrays),  the  probability  of  at  least  K  detections  out 
of  M  is  given  by  the  summed  binomial  distribution 


l*(>K) 


M 


l  (j)p‘(l-p) 
i  *K 


M-i 


(1) 


This  formula  is  used  to  evaluate  either  the  falsc-alarm 
rate  or  the  signal  detection  probability,  and  p  is 
accordingly  computed  from  a  central  l:  or  a  non-central 
P,  respectively.  Subroutines  to  do  these  calculations 
arc  given  in  the  Appendix.  Another  assumption  we  have 
made  in  our  calculations  here  is  that  the  frequency 
window  is  1/2  hertz,  so  that  Nj  ■  T,  i.c.  the  number  of 
degrees  of  freedom  of  the  numerator  equals  the  time 
window  in  seconds,  an  approximation  to  a  system  already 
in  operation  (Blandford,  1970).  In  any  case,  the  non¬ 
centrality  parameter  is  A  ■  N  Nj  (S/N)^,  where  (S/N)  is 
the  signal/noisc  ratio  on  each  clement  and  N  is  the 
number  of  elements/subarray.  In  expressing  our  falsc- 
alarm  rate  in  terms  of  falsc-alarms/day  we  have  assumed 
86400/T  samples  per  day,  although  if  an  overlapping  or 
sliding  window  is  used  (the  usual  case),  there  will 
presumably  be  an  additional  as-yct-undetermined  factor. 

In  view  of  the  steepness  of  the  operating  curves 
obtained,  this  is  probably  not  a  crucial  point. 

We  have  evaluated  this  detector  for  certain  con¬ 
figurations  that  might  be  used  at  N’ORSAK:  6  elements/ 
subarray  (Nj  ■  3  and  *  15)  and  7,  13,  and  22  subarrays. 
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Figures  6  through  8  show  the  effect  of  choice  of  K  on 
signal  detection  vs  false-alarn  rate  for  a  fixed  signal 
strength.  Figures  9  through  11  arc  operating  curves  for 
the  optimum  K  detectors  for  the  same  configurations. 
Table  II  summarizes  these  results  in  terms  of  signal/ 
noise  ratios  on  subarray  beams  for  90  per  cent  detec¬ 
tion  probability  with  one  false-alarra/day  and  compares 
them  with  the  performance  of  a  single  subarray  and  with 
the  ideal  performance  of  a  single  F  detector  operating 
on  the  whole  array  ("total  beam"),  given  in  Table  I. 
Also  included  in  Tublcs  I  and  II  are  similar  data  for 
groups  of  12-element  subarrays  (Nj  »  3  and  Nj  ■  33). 
figures  4  and  5  give  operating  curves  for  single  sub¬ 
arrays.  (Note  that  element  signal/noisc  ratios  arc 
derived  from  beam  signal/noisc  ratios  by  dividing  by 
the  square  root  of  the  number  of  elements  in  the  sub¬ 
array.)  Table  II  shows  that  the  optimum  K  detector  is 
significantly  superior  to  a  single  subarray,  and  in 
the  ease  where  the  total  array  is  so  large  that  sig¬ 
nificant  loss  of  signal  coherence  occurs,  it  may  well 
perform  better  than  a  single  F  detector  operating  on 
the  total  beam. 

In  view  of  the  result  obtained  above,  that  using 
a  time  window  somewhat  longer  than  the  three  seconds 
average  signal  length  gave  better  results  with  a  single 
detector,  we  tried  using  a  longer  time  window  on  the 
composite  detector.  The  results,  shown  in  Table  III, 
are  negative,  even  in  the  ease  of  only  three  elements/ 
subarray.  Apparently  the  composite  detector  is  so  much 
mere  sensitive  to  signal/noise  ratio  than  to  degrees  of 
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freedom  that  the  improvement  in  noise  estimate  obtained 
with  a  longer  window  is  insufficient  to  offset  the 
degradation  in  effective  signal/noise  ratio.  We  have 
assumed  a  signal  exactly  three  seconds  long,  so  the 
effective  S/N  ratio  will  be  degraded  by  a  factor  /3/T. 

Of  course  in  practice  there  will  usually  be  some  contri¬ 
bution  to  the  signal  after  three  seconds,  but  it  is 
probably  still  true  that  the  best  time  window  for  a 
network  of  F  detectors  corresponds  to  the  effective 
signal  duration,  even  though  the  best  time  window  for 
an  individual  F  detector  may  be  somewhat  longer.  Thus 
in  future  systems  there  may  have  to  be  a  compromise 
between  the  role  of  a  station  as  a  station  and  its  role 
as  part  of  a  network,  but  the  difference,  in  any  case, 
is  not  very  great. 

For  network  applications  there  is  another  consid¬ 
eration  that  has  a  possible  bearing  on  the  choice  of 
time  window  and  on  the  way  decisions  are  reached.  A 
calculation  based  on  (1)  shows  that  for  a  22-station 
system,  a  network  falsc-alarm  rate  of  1  per  day 
implies  an  average  falsc-alarm  rate  of  ~4600/day  for 
each  individual  station.  Furthermore,  decreasing  the 
network  false-alarm  rate  to  10"*Vday  decreases  the 
station  false-alarm  rate  only  slightly,  to  ~1300/day, 
due  to  the  K'th  power  in  (1).  Thus  the  individual 
station  detections,  b£  themselves .  would  be  relatively 
useless  for  decision  making,  at  the  thresholds  used 
for  network  detection.  Also,  experimentally  determining 
the  proper  delays  to  correlate  individual  "detections" 
would  seem  impractical  because  of  the  large  number  of 


-9- 


combinations.  The  best  approach  would  seem  to  be  a 
fully  automatic  detector  in  which  "all  possible  beams" 
are  scanned  to  cover  a  given  area.  In  practice  this 
would  merely  involve  scanning  a  table  of  delay  times 
and  counting  votes  over  sections  of  data  corresponding 
to  the  greatest  and  least  delays  for  each  station  for 
that  region.  It  does  not  seem  at  all  unlikely  that  a 
region  the  size  of  Russia  could  be  monitored  in  real 
time  with  existing  computers.  With  an  efficient  pro¬ 
cedure  mucii  larger  areas  could  probably  be  scanned.  The 
size  of  each  beam  and  hence  the  number  necessary  to 
cover  a  given  area,  would  be  influenced  by  the  time 
window  used:  a  longer  time  window  would  reduce  the 
number  required.  Notice  that  automatic  network  detectors 
necessarily  give  approximate  locations  simultaneously 
with  the  detections.  Most  of  these  remarks  also  apply 
to  the  multi-array  F  detector  described  below. 

We  should  also  point  out  that  false-alarm  rates 
quoted  in  this  report  are  per  beam  and  must  be  multi¬ 
plied  by  the  number  of  beams  to  give  the  total  network 
false-alarm  rate.  For  a  system  operating  a  thousand 
beams,  the  rates  quoted  here  are  certainly  too  high. 
False-alarm  rates  down  to  10  /day  can  be  read  directly 
from  the  operating  curve;,  and  linear  extrapolation  is 
satisfactory.  For  a  22  station  system,  threshold  mag¬ 
nitude  increases  ~ , 6dB/decade  decrease  in  false-alarm 
rate.  This  represents  another  way  in  which  time  window 
length  can  influence  detection  thresholds  in  real 
systems . 

A  variation  of  the  composite  detector  was  also 
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investigated:  on  receipt  of  at  least  k^  subarray  detec¬ 
tions  at  threshold  Fp  the  threshold  on  the  remaining 
subarrays  is  lowered  to  F2  and  a  total  of  at  least  k2 
detections  is  then  required  in  order  to  declare  a 
detection  by  the  system.  Writing  P(jJp)  to  represent 
the  probability  of  exactly  k  events  out  of  N  when  indi¬ 
vidual  events  have  probability  p,  a  shorthand  for  the 
binominal  distribution,  the  probability  of  at  least  k2 
subarrays  detecting  by  this  scheme  is 


P(>k2) 


P( 


N 

— k2 


Pi) 


V1 


P(ilpl')  P<:>(k2-i)  lp2“pl) 


i=k 


(2) 


where  p^  is  the  subarray  detection  probability  at  thres¬ 
hold  Fp  p2  is  the  probability  at  threshold  Fp  and 
p2-p^  is  the  fraction  lying  between  F2  and  Fp  (We  use 
the  symbol  "  >  "  to  mean  "at  least".)  Note  that 

—  P*->k  I P i )  •  This  scheme  was  always  worse  than 
the  previous,  and  much  simpler,  detector.  It  gave  com¬ 
parable  results  only  in  the  limits  k^  -*■  k2  and/or 
P1  p2*  w*i:'-ch  reduce  to  the  previous  case.  (False- 
alarm  rate  was  always  held  constant;  the  ratio  F2/F^ 
was  specified  and  F^  determined  by  an  inversion;)  The 
reason  for  the  poorer  performance  is  probably  that  the 
effective  number  of  subarrays  is  reduced  by  using  a 
lower  K  value  and  by  exaggerating  the  weak  subarrays: 
the  statistical  "inertia"  of  the  system  is  reduced  by 
the  manipulations.  Schemes  of  this  type  appear  to  be 
extremely  unpromising  and  the  idea  has  not  been  pursued 
further.  It  might  still  be  worthwhile  to  review 
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non-detecting  stations  for  arrival  time  information  in 
order  to  obtain  a  better  least-squares  location,  although 
it  may  very  well  be  that  using  the  arrival  times  on  such 
stations  would  not  actually  produce  a  better  location. 
This  is  particularly  apt  to  be  true  for  large  K,  where 
one  is  practically  guaranteed  a  good  location  for  all 
detections  anyway. 
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A  MULTI -ARRAY  F  DETECTOR 


There  is  another  approach  possible  when  signals  are 
not  the  same  on  all  subarrays,  which  is  a  special  case 
of  a  technique  of  Shumway  (1970).  We  consider  the  model 
representing  the  output  from  the  i'th  sensor  on  the 
j  ' th  subarray  as 


yn(t) 


aisi(t) 


where  j  =  1, 


M;  1  =  1, 


•  •  •  9 


N;  and  t  =  0,  1, 


L-l. 


to  be  noimal,  stationary,  and 


We  assume  the  noise  n.^0 
uncorrelated  between  sensors,  and  the  signals  Sj  to  be 
the  same  on  all  sensors  of  a  given  subarray.  We  further 
assume  the  Sj  all  to  have  the  same  rms  value  over  time 
and  explicitly  allow  for  different  size  signals  by 
inclusion  of  the  factors  a^ ,  which  we  assume  are  known 
a  priori.  Then,  following  Shumway  (1970),  the  maximum- 
likelihood  estimate  for  the  signal  on  the  j'th  subarray 
is  just 


V°  ’  aj:1  yj-U) 

where  the  dot  signifies  the  subscript  over  which  the 
mean  is  taken. 

In  order  to  test  the  hypothesis  that  there  is  no 
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signal  present  on  any  subarray,  i.e.,  that  all  S.(t)  »  0, 
we  could  transform  to  the  frequency  domain  and  form  the 
ratio  of  beam  power  to  residual  power  as  our  test  sta¬ 
tistic.  But  if  t lie  data  arc  band-limited,  to  bandwidth 
B,  over  which  the  noise  spectrum  is  constant,  wc  can 
instead  use  an  approximation  to  the  F  statistic  given 
by 


1 2BTM, 2BTM(N- 1)  = 


(n-i)n  IjItyj5co 
lj  lt  II,  yj-(t)  -  n  yjfct) j 


(3) 


where  B  is  the  filter  bcadwidth  in  hertz,  T  the  time 
window  in  seconds,  M  the  number  of  subarrays,  and  N 
the  number  of  eloments/subarray.  In  the  presence  of 
signal,  the  non-centrality  parameter  is 


x<u0) 


2BTN  I  a.2|S(u  )| 
j  "1  ' 


N2(u,0) 


(4) 


where 


,  ij.t  *jz(t) 


l^(w0)  I  ” 


BT  y.  a. 

L)  J 
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is  the  approximate  mean  signal  power  within  the  frequency 
band,  and  N  (<jq)  is  the  noise  power.  We  include  the  fil¬ 
ter  center  frequency  uQ  in  (4)  in  order  to  emphasize  the 
dependence  on  filter  characteristics;  in  all  of  the  work 
in  this  report,  we  assume  a  filter  similar  to  that  used 
at  Tl:0,  as  described  in  Blandford  (1970). 

W'c  can  abbreviate  (4)  to 


X  -  N  N.  (  1  a . 2/M)  (S/N)2 

1  j-1  J 


(5) 


where  Nj  =  2BTM  is  the  number  of  degrees  of  freedom  of 
the  numerator  and  (S/N)  is  the  normalized  signal/noisc 
ratio,  defined,  when  aj  ■  1,  as  in  Blandford  (1970). 

When  all  aj  ■  1,  the  expression  for  X  has  the  same 
form  as  previously,  although  the  statistic  is,  in 
general,  different.  When  M  ■  1,  the  detector  reduces 
to  the  single-array  F  detector. 

In  order  to  evaluate  the  performance  of  this  detec¬ 
tor  on  an  array  such  as  NORSAR,  where  the  signals  arc 
generally  of  the  same  size  on  all  subarrays,  we  set 
all  aj  ■  1.  Table  IV  gives  the  results  for  some  of  the 
same  configurations  for  which  the  previous  detector 
was  evaluated.  The  multi-array  F  detector  does  somewhat 
better  but  is  still  worse  than  the  ideal  total  array 
beam.  Figures  12  through  14  give  operating  curves  for 
this  detector.  Results  obtained  for  a  six-second  time 
window  were  worse,  indicating  that  there,  too,  the 
optimum  length  is  about  three  seconds  for  a  three 
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second  signal.  Tabic  V  gives  operating  thresholds  (I 
values)  for  all  the  systems  described. 
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The  multi-array  1  detector,  as  presented  here, 
could  also  he  used  on  a  worldwide  network  of  stations. 

We  have  so  far  assumed  the  noise,  as  well  as  the  signal, 
to  have  the  same  rms  amplitude  on  all  arrays.  When  this 
is  net  the  case,  we  first  normalise  the  traces  hy  a 
long-term  average  of  the  noise  estimates.  The  average 
should  he  long  enough  not  to  affect  the  statistics  and 
short  enough  to  follow  daily  variations  in  the  noise 
level.  In  practice,  a  running  average  of  the  estimated 
residual  noise  power  from  the  preceding  hour  or  so 
would  probably  be  satisfactory.  The  theory  assumes  the 
noise  to  be  identically  distributed  on  all  channels, 
but  amplitude  variations  can  be  treated  as  an  amplifier 
gain  error.  Notice  that  in  forming  the  network  test 
statistic  (3),  the  sums  over  i  and  t  can  be  done  at 
each  station  prior  to  transmission  to  the  central  sta¬ 
tion,  so  that  the  only  quantities  that  need  to  be 
transmitted  arc  the  contributions  to  the  j  sums  in  the 
numerator  and  denominator,  i.c.,  just  two  numbers.  (The 
contributions  to  the  denominator  will  average  close  to 
1,  due  to  the  noise  normal i ration,  but  the  fluctuations 
arc  important.)  The  central  station  then  has  to  sum 
contributions  from  windows  with  the  appropriate  time 
delays,  divide,  and  compare  with  the  threshold. 

For  purposes  of  analyzing  the  performance  of  such 
a  network,  the  amplitude  factors  3j  would  be  given  by 
the  distance-amplitude  corrections  for  each  station 
for  a  given  epicenter  location,  divided  by  the  rms  noise. 
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Relation  (5)  says  that  X,  and  hence  the  detection  capa¬ 
bility,  depends  on  the  root-mean-squaro  of  the  amplitude 
factors.  Thus  curves  such  as  given  in  figures  12  through 
can  be  used  here  also  by  multiplying  the  signal/noise 
ratio  indicated  in  the  Figures  by  a^(J  a ( “/M) " ‘  ^  to 
obtain  the  signal/noise  ratio  required  on  each  station. 
Note  that  the  false-alarra  rate  does  not  depend  on  the 
a  ,  so  Table  V  can  also  be  used  for  this  case. 

The  composite  (>K)  detector  can  also  be  used  for 
a  worldwide  network,  and  since  there  are  no  special 
requirements  on  the  noise  at  different  stations,  no 
normalization  need  be  done.  In  addition,  the  summation 
at  the  central  station  should  he  somewhat  faster  than 
in  the  case  of  the  multi-array  F  detector,  so  that 
more  beams  could  be  formed  in  real  time.  Also,  the  only 
information  that  needs  to  be  transmitted  to  the  central 
station  is  a  1  or  a  0.  The  difficulty  with  this  detector 
lies  in  the  analysis  of  its  performance,  which  is  trouble¬ 
some  but  not  insurmountable.  The  first  thing  to  notice  is 
that,  with  unequal  signal  strengths,  the  probabilities 
of  detection  on  each  individual  station,  pj ,  will  no 
longer  be  equal  as  required  by  (1).  Thus  to  compute 
l’(>K),  a  more  complicated  procedure  must  be  used  (Wirth, 
1971).  The  falsc-nlarn  rate  can  be  computed  as  before 
if  all  the  stations  have  the  same  number  of  elements, 
although  this  is  not  a  requirement  for  the  detector  to 
*orh.  (One  could  require  that  each  station  have  the 
same  fa lso-alarm  rate,  and  then  set  different  thresholds 
1-v  for  stations  with  different  numbers  of  elementSi  or 
degrees  of  freedom.  Numerically  this  wouid  b«  easy  to 
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do  und  might  also  possess  some  statistical  idvantagc 
over  the  scheme  of  using  the  same  I  threshold  on  each 
station,  although  tins  has  not  been  investigated.)  The 
principal  difficulty  lies  in  the  fact  that  the  optimum 
choice  of  k  depends  on  the  relative  signal  strengths  at 
the  different  stations,  as  will  he  shown  hclow.  Thus 
for  optimum  performance,  different  choices  of  K  will  in 
general  he  required  for  monitoring  different  regions  of 
the  earth.  There  would  he  no  particular  problem  in 
implementing  this,  since  k  is  just  the  minimum  number  of 
individual  "detections"  required  in  order  for  the  cen¬ 
tral  collecting  station  to  proclaim  a  network  detection, 
and  different  numbers  could  easily  he  programmed  for 
monitoring  different  regions  of  the  earth.  A  complica¬ 
tion  arises  because  t lie  falsc-alarm  rate  also  depends  on 
k,  so  one  would  cither  have  to  allow  the  falsc-alarm 
rate  to  vary  for  different  epicenter  regions  or  else 
specify  different  station  thresholds  also.  Again,  this 
would  not  be  difficult  to  implement.  Thus  the  composite 
(>k)  detector  is  somewhat  more  flexible,  slightly  less 
sensitive,  and  requires  a  lot  more  individual  tailoring 
than  the  multi-array  F:  detector  for  network  applications. 

To  investigate  the  behavior  of  the  optimum  k  for 
differing  signal  levels,  we  consider  a  simple  case  in 
which  there  arc  only  two  different  signal  levels.  We 
take  Mj  stations  with  beam  signal/noisc  ratio  (S/N) j 
and  stations  with  ratio  (S/N),,  such  that  M  »  Mj  ♦  M? 
is  constant.  Overall  falsc-alarm  rate  is  held  constant, 
figure  15  is  a  contour  plot  cf  the  optimun  k  values, 
holding  (S/N) 2  constant  and  allowing  (S/N) j  to  vary 
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geometrically.  (The  surface  should  be  considered  to  be 
composed  of  terraces.  Sec  Kirth,  1971,  for  a  description 
of  the  contour  algorithm  used.)  The  figure  shows  that 
the  general  effect  of  having  a  few  stations  with  large 
signaj./noise  ratio  is  to  decrease  the  value  of  K.  The 
figure  also  makes  it  obvious  that  guessing  the  optimum 
k  would  be  a  little  difficult,  even  in  this  simpliicied 
case.  Given  an  actual  network,  however,  it  would  not  be 
difficult  to  calculate  the  optimum  strategies.  Figure 
»6  gives  contours  of  the  corresponding  probability  of 
detection  for  the  optimum  choices  of  K.  The  increase 
towards  the  upper  right  corner  is  not  really  surprising. 

An  interesting  comparison  can  be  made  with  the 
multi-array  F  detector  by  holding  the  rms  of  all  the 
signals  constant,  instead  of  (S/N)^.  Since  the  F 
detector  is  sensitive  only  to  the  root-mcan-squarc  of 
the  signal/noise  ratios  on  all  the  stations,  the  prob¬ 
ability  of  detection  should  be  constant.  Figures  17  and 
18  arc  the  corresponding  plots  for  this  constraint.  The 
rms  signal/noisc  ratio  is  1.224,  for  which  the  multi¬ 
array  F  detector  has  a  90  per  cent  probability  of 
detection.  The  behavior  of  K  is  very  similar  to  the 
preceding  case,  indicating  that  it  is  sensitive  primarily 
to  the  ratios  of  signal  levels.  The  behavior  of  the 
probability  of  detection  for  the  optimum  choice  of  K 
(Figure  18)  is,  however,  not  entirely  expected.  It 
should  be  remembered  that  the  multi-array  F  detector 
should  have  90  per  cent  probability  of  detection  over 
the  whole  plot.  This  points  up  another  difference 
between  the  two  detectors:  while  for  equal  signal 
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levels  the  composite  (>K)  detector  is  nearly  as  good  as 
the  multi-array  !•*  detector,  for  disparate  signals  it 
docs  increasingly  worse.  In  particular,  figure  18  shows 
that  the  presence  of  a  few  stations  with  large  signal/ 
noise  ratio  is  especially  damaging  in  comparison.  This 
should  be  kept  in  mind  when  making  a  choice  between  the 
two  detectors  for  any  particular  system,  however  the 
difference  may  not  be  as  great  as  it  seems.  The  66  per 
cent  contour  represents  a  1  dB  difference  in  signal 
detection  thresholds  between  the  two  detectors  (Table  IV) 
and  the  30  per  cent  contour  probably  represents  only  an 
additional  1  dB,  judging  from  the  operating  curves. 

We  have  so  far  based  our  analysis  on  the  assumption 
that  all  stations  were  to  be  retained,  no  matter  how 
poor.  This  is  not  necessarily  the  way  things  arc  done 
in  real  life,  and  a  little  reflection  shows  that  this 
is  not  always  wise,  for  the  composite  detector  it  is 
easy  to  see  why:  very  poor  stations  contribute  only 
false  alarms  and  thus  represent  a  ncgativ'  asset.  The 
same  thing  is  true  of  the  multi-array  F  detector,  but 
the  reason  is  not  quite  so  obvious.  The  work  on  optimum 
time  windows  for  this  detector  shows  that  increasing 
the  degrees  of  freedom  while  holding  X  constant  results 
in  poorer  performance  unless  and  arc  very  small. 
Equations  (3)  and  (4)  say  that  this  is  exactly  what 
happens  when  one  adds  very  weak  stations  to  the  network 

(aM+l<<c  l  aj2)* 

The  computations  done  above  for  the  composite 
detector  were  repeated,  throwing  out  the  poor  stations. 
Figures  19  and  20  are  analogous  to  Figures  15  and  16 
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and  give  the  contours  of  optimum  K  and  the  probabilities 
of  detection  for  varying  numbers  of  strong  stations  (M) 
and  signal/noise  ratio  (S/N).  Calculations  were  also  done 
extending  the  graphs  downwards,  i.e.,  for  the  case  of  a 
few  weak  stations.  These  calculations  show  that  it  is 
generally  better  to  discard  stations  whose  signal/noise 
ratio  is  less  than  about  .4  -  .6  times  the  signal/noise 
ratio  on  the  strong  stations,  if  M  >  3.  (For  example, 
comparing  Figures  16  and  20  along  the  lines  or  M  =  6 
shows  a  crossover  near  (S/N) (S/N) ^  ~  2.0,  above  which 
it  is  better  to  drop  the  weak  stations.)  Considerations 
about  the  minimum  desirable  number  of  stations  in  the 
network  for  determining  locations  will  also  play  a  part 
here.  If  half  the  stations  (11)  have  only  1/4  the 
signal/noise  ratio  on  the  other  half  (which  is  1.5), 
then  the  probability  of  detection  is  -60  per  cent  if 
the  weak  stations  are  dropped,  compared  with  -30  per 
cent  if  they  are  retained,  which  corresponds  to  a 
"threshold  magnitude"  change  of  about  1  dB.  (In  Figures 
16  and  20,  the  greatest  difference  would  be  in  the  upper 
right  corner,  but  since  the  probabilities  are  close  to 
1  here  anyway,  the  difference  is  not  as  great  as  in 
the  example  cited.)  Quantitative  studies  of  the  multi¬ 
array  F  detector  for  this  case  are  being  carried  out  in 
the  context  of  determining  the  optimum  filter,  which  is 
an  analogous  problem,  as  pointed  out  by  Blandford,  These 
will  be  reported  in  the  future.  Preliminary  results 
indicate  a  rejection  criterion  similar  to  that  for  the 
composite  detector. 
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CONCLUSIONS 


1.  A  signal  time  window  longer  than  3  seconds 
improves  the  detection  capability  of  F  detectors  on 
small  arrays,  resulting  in  performance  closely  approach¬ 
ing  that  of  a  power  detector,  without  degrading  the 
desirable  features  of  the  F  detector. 

2.  These  longer  windows  are  a  disadvantage,  however, 
when  the  arrays  are  used  in  a  detection  network.  Quali¬ 
tatively  this  is  because  more  than  one  station  must 
trigger  for  a  detection,  so  that  a  false  alarm  on  only 
one  due  to  statistical  fluctuation  is  not  serious.  In 
practice  both  the  short  and  long  windows  might  be  used, 
the  long  for  detection  at  each  station  independently; 

and  the  short  for  the  final  decision  in  combining  the 
station  detections.  After  detection,  a  resurvey  of  the 
data  could  be  performed  for  location  purposes. 

3.  A  full  22-subarray  NORSAR  can  perform  within 
5.4  dB  of  its  theoretical  N  '  performance  for  perfect 
signal  correlation  by  using  independent  F  detectors  on 
each  6-element  subarray  and  declaring  a  detection  with 
12  (the  optimum  number)  or  more  subarray  detections. 

A  multi-array  F  detector  would  lose  only  4.3  dB.  An 
iterative  technique  in  which  a  first  detection  is  made 
with  a  higher  threshold  on  fewer  than  K  subarrays, 
followed  by  detection  using  a  lower  threshold  on  more 
than  K,  is  uniformly  worse  than  the  simple  method. 

4.  If  some  stations  have  higher  signal-to-noise 
values  than  others,  the  advantage  of  the  multi -array 
detector  over  the  composite  detector  becomes  somewhat 
greater . 
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5.  The  multi-array  detector  requires  slightly  more 
data  transmission,  but  the  decision  strategy  for  the 
composite  detector  is  more  difficult  to  calculate. 

6.  The  different  network  detectors  can  be  imple¬ 
mented  for  a  large-aperture  array,  or  a  real  world-wide 
network,  by  straightforward  extrapolation  of  the  tech¬ 
niques  presented  in  this  paper. 
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TABLIi  I 


Detection  levels  for  total  beams,  assuming  signals  per¬ 
fectly  correlated.  (Signal/noisc  ratios  required  for  90* 
detection  at  1  false-alarm/day ,  with  3  second  time  window 
and  Nj  *  3). 
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WINDOW  IN  SECONDS 

Figure  1.  Signal- to-noise  values  required  at  the  seismometer 
level  for  a  false  alarm  rate  of  0.1  per  day  and  probability 
of  detection  of  0.9  as  a  function  of  the  signal  window  in 
seconds.  Curves  are  shown  for  arrays  ranging  from  2  to  31 
elements,  and  another  curve  lias  been  drawn  through  the  minima 


FALSE  ALARMS/DAY 


PROBABILITY  OF  DETECTION  -  % 


Figure  4,  Probability  of  detection  versus  beam  signal/noise 
single  subarray  of  6  elements. 
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Figure  6.  Probability  of  at  least  K  subarrays  out  of  7  detecting 
(6  elements/subarray) . 
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Figure  7.  Probability  of  at  least  K  subarrays  out  of  13 
detecting  (6  elements/subarray). 
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Figure  9.  Probability  of  at  least  5  subarrays  out  of  7  detecting 
versus  beam  signal/noise  (K  optimum,  6  elements/subarray) . 
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Figure  10.  Probability  of  at  least  8  subarrays  out  of  13 
detecting  versus  beam  signal/noise  (K  optimum,  6  elements/ 
subarray) . 
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Figure  11.  Probability  of  at  least  12  subarrays  out  of  22 
detecting  versus  beam  signal/noise  (K  optimum,  6  elements/ 
subarray) . 
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Figure  12.  Probability 
versus  subarray  beam  si 
sub  array) . 
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Figure  13.  Probability  of  detection,  multi-array  F  detector, 
versus  subarray  beam  signal/noise  (13  subarrays,  6  element'./ 
subarray) . 
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STATISTICAL  CALCULATION'S 


(if)  1 


NONCl.NTKAl.  I;:  Subroutines  have  been  written  to  evaluate 
the  cumulative  distribution  Function  of  the  noiiceutrai  I 
distribution  for  all  parities  of  Vj  and  y,.  The  necessary 
formulas  have  been  taken  from  Abramowitr  and  St  cguu  (I960), 
with  corrections,  i lie t r  formula  .0.0. IB,  k i v i ii^t  the  rela¬ 
tionship  between  the  central  and  non-central  distributions, 
is  incorrect,  as  pointed  out  by  Itlnndford  (1970).  Ibe 
correct  form  is 


<1) 


The  upper  integral  Q  proh(u  ^  I  )  is  the  one  evaluated 
here,  and  the  relationship  is  Q  ■  1  -  1*.  It  is  easy  to 
sec  that  (1)  still  holds  if  the  l'*s  on  both  sides  are 
replaced  by  Q's.  This  formula  is  used  for  all  three  cases. 
The  sum  always  converges  for  finite  X,  however,  overflows 
may  occur  If  e'^  exceeds  the  range  of  floating-point  vari¬ 
ables.  On  the  CDC  1604  computer,  this  restricts  the  non- 
centrality  parameter  to  1  <  1416.  The  mux. mum  of  2000  terms 
of  (1)  taken  in  the  subroutines  should  he  adequate  within 
this  range  of  X  for  an  accuracy  of  at  least  six  significant 
figures  in  Q. 

Observe  that  the  third  parameter  (v,)  of  Q  on  the  RIIS 

•» 

of  (1)  docs  not  depend  on  j.  Use  is  made  of  this  fact  in 
choosing  formulas  for  QCFlvpV^)  in  order  to  minimize 
computations.  It  has  been  possible  in  two  cases  to  obtain 
forms  which  require  only  ~M  operations,  where  M  is  the 
number  of  terms  taken  in  the  infinite  scries,  rather  than 


-Al- 


M"  as  might  have  teen  expected.  In  the  third  case,  v, 
even,  the  infinite  sum  can  he  done  explicitly.  These 
cconemies  are  possible  because  the  parameter  x  in  all 
the  expansions  for  Q,  defined  by 


does  not  depend  on  j. 

v j  even:  For  this  case,  formula  26.6.4  is  used  for  Q: 

v,/2  v?  v,(v2*2)  2 

QCFlVpVj)  *  x  *  (1*7^  (1-x)  ♦  -£7^ -  (l-x)z  *  ••• 

(3) 


,Cyi)—  (V;>V|-t)  *,«*•. 

^ • 4 • • • f j  t J 


where  x  is  defined  by  (2).  Since  the  only  dependence  on 
Vj  (and  hence  on  j)  is  in  the  number  of  terms,  each 
additional  term  in  the  infinite  series  in  (1)  requires 
only  one  additional  term  in  (3),  supporting  the  claim 
made  in  the  paragraph  above.  Listings  of  all  the  sub¬ 
routines  (QNCF)  arc  given  at  the  end  of  this  section. 

v j  and  v2  odd:  Formula  26.6.8  can,  with  some  manipu¬ 
lation,  be  written  in  the  more  convenient  form 


(in  3 
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Q0:IVv2> 


■ 


♦  •  •  • 


-  U  -  x  cot  $  1 1  ♦  i  x 


2  •  4  •  •  •  ( 3)  (w2*3)/2.j 
*  3»5*»»r,v,-2)  X 


2*4***(v,-l)  (v,*1)/2 

4  *•.  *  x  6  cot 

n  (Vj-ZJ  * 


v,*l 

(1  *  -4—  (1-x) 


(U 


•  •  •  • 


(v2*l) (v2*J)*»»(vl*v2-4)  t  ^(Vj-J)/2, 
3*  5* •• CVj-2)  ) 


where 


<  i  a  retail  /v2/Vj^ 

and  x  is  defined  hy  (2).  Note  that,  once  again,  the  only 
dependence  on  Vj  is  in  the  number  of  terms  (of  the  second 
scries),  so  that  an  efficient  calculation  is  possible  here 
also.  Use  is  also  aadc  of  the  fact  that  the  first  series 
can  be  rcaovcd  froa  the  iminitc  sua  since  it  is  indepen¬ 
dent  of  j. 

v2  even:  formula  26.6.22  is,  unfortunately,  incorrect, 
nor  did  several  obvious  modifications  of  the  formula 
yield  correct  answers.  A  similar  formula  was  derived  by 
Wirth  (1971),  starting  from  26.6,5  written  in  the  form 


P(F|vj,v2) 


Vj/2 

(1-x)  1  (1 


V1 

T~ 


x  ♦ 


(vi*2) 

— rn— 


x2  ♦ 


.  .  . 


♦ 


V1^'V2)***  (vi  +  v2‘4) 
2*4» • • (v2*2) 


v2/2-1 
x  2  J 


C>) 
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The  derivation  will  not  be  repeated  hero,  but  involves 
putting  (5)  into  (1),  interchanging  the  order  of  summa¬ 
tions,  applying  Rummer's  transform  to  the  confluent 
hypergeometric  series,  and  deriving  a  recursion  for  the 
resulting  finite  series.  The  result  is  a  closed  expres¬ 
sion  for  Q  which  is  very  convenient  for  computation 

v  /i 

Q(P|v, ,»,,»)  •  i  -  c*>x/2  (1-x)  1  £  T, 

1  •  l-o  * 


where 


T1  *  y  Ivj  ♦  X(l-x)) 

Ti  -  7T  Uv1*4i-4*X(1-X)JT4.1-  (v^2i-4)xTiw2) 


and  x  is  defined  by  (2). 

The  three  subroutines  were  checked  against  tables 
of  the  central  F  distribution  in  Abramowitz  and  Stegun 
(1966)  with  \  ■  0.  The  routines  for  Vj  even  and  V2  even 
were  checked  against  each  other  and  against  routines 
which  had  been  written  independently  by  Blandford  (1970) 
from  formulas  26.6.6  and  26.6.7.  In  addition,  all 
routines  were  checked  against  the  Pearson  and  Hartley 
charts  reproduced  in  Schcffc  (1959),  relating 


noncontrulity  parameters  according  to  A  «  The 

agreement  in  ail  wascs  was  within  the  limits  of  road¬ 
ability  or  table  accuracy,  or  the  normal  limits  of 
single-precision  computation.  The  nominal  accuracy  of 
the  first  two  routines  is  about  six  digits,  due  to  the 
truncation  of  the  infinite  series. 

INVURSli  F:  A  very  useful  routine  has  also  been  developed 
for  computing  the  inverse  of  the  central  F  distribution 
function,  F(Q»Vj,v2),  which  gives  the  threshold  level 
corresponding  to  a  given  false-alarm  rate.  This  routine 
has  made  possible  the  automatic  plotting  of  detection 
probability  directly  in  terms  of  false-alarm  rate. 
Newton's  iteration 


F 


i*l 


Q(Pj)  •  Q 

Q'NV 


is  used,  with  the  previously  described  routines  (QNCF) 
being  used  to  compute  Q(F).  The  derivative  is  obtained 
from  26.6.1  by  differentiation 


Q'(F) 


V2 


(rF) _ 

FBOiVj  ,*|V2)  (l*rF)  <VV2)/2 


where  r  i  Vj/v2  and  B  is  the  beta  function: 

B(a,b)  =  r(a)r(b)/r(a»b).  Starting  values  are  derived 
from  26.6.16  and  26.5.22  over  part  of  the  range  of 


arguments.  The  inverse  of  the  normal  cumulant  required 
for  this  formula  is  provided  by  a  routine  (QUANTF) 
derived  from  a  Hastings  expansion,  formula  26.2.23.  Over 
the  rest  of  the  range,  a  starting  value  is  obtained  from 
the  asymptotic  form 


Q(f|v1,v2) 


2 

v2B  (l  +  rF)V2/2 


for  rF  >>  1 


from  which 


F 


1 

r 


[( 


2/v. 


i] 


where  arguments  of  B  are  omitted  for  simplicity.  (This 
formula  is  exact  for  ■  2.)  The  resulting  routine 
(FINV)  converges  quite  well  over  a  broad  range  of  input 
parameters.  In  less  than  10  iterations  it  will  match 
the  input  value  of  Q  to  at  least  5  digits  relative 
accuracy.  (Note,  however,  that  the  second  two  versions 
of  QNCF  have  a  lower  limit  on  absolute  precision  of 
~10"10  due  to  rounding.  Thus  the  relative  accuracy  will 
be  reduced  for  Q  <<  10"^.) 

An  inverse  of  the  noncentral  F  might  also  be  useful, 
X  (F,Q,v^,y>2)  ,  giving  the  noncentrality  parameter  in 
terms  of  the  probability,  i.e.  the  signal/noise  ratio 
required  for  a  given  probability  of  detection.  The 
derivative  for  a  Newton's  iteration  (7)  is  easy  to 
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obtain  from  (1)  : 


3QCF| v1,v2,A) 

- 5T - 


Q(F|vlfv2,»)] 


but  it  is  not  so  easy  to  see  how  to  derive  a  starting 
value.  The  approximations  in  Abramow’  z  and  Stegun  are 
impractical  for  this  purpose.  Actually,  such  a  routine 
has  not  been  needed  here.  Signal/noise  thresholds  quoted 
in  the  text  were  obtained  graphically  by  linear  inter¬ 
polation  with  respect  to  (S/N).  The  high  linearity  of 
the  curves  with  respect  to  both  (S/N)  and  false-alarm 
rate  may  furnish, a  clue  to  obtaining  a  starting  value 
for  the  automatic  calculation,  should  this  approach 
seem  desirable  in  the  future. 

BINOMIAL:  The  routine  for  computing  the  summed  binomial 
distribution,  equation  (1)  of  the  text,  is  completely 
straightforward.  A  listing  of  the  subroutine  (POFK)  is 
given  below. 

INVERSE  BINOMIAL:  The  inverse  of  the  binomial  distri¬ 
bution  is  used  with  the  inverse  of  the  F  distribution 
to  obtain  the  threshold  in  terms  of  the  false-alarm 
rate  for  the  composite  (>K)  detector.  A  Newton's 
iteration  (7)  is  used  with  a  modified  derivative.  The 
derivative  of  P(>k)  is  easily  obtained  from  (1)  of  the 
text 
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0  f)  8 


(8) 


dP(>k) 

Up 


=  k  (“) 


P^U-p)^ 


and  it  is  easily  computed,  however,  the  iteration  did 
not  converge  as  well  as  might  have  been  expected  when 
this  derivative  was  used  with  a  starting  value  derived 
from  an  approximation  to  P(>k).  It  is  easily  shown  that 
the  lowest  term  of  P (>k)  is  the  largest  if  p  <  (k+l)/(M+l), 
so  that 


P(>k)  s  (jj)  pk(l-p)M_k 


from  which  a  starting  value  p  »  [P  (>k)  /  (^)  ]  ^k  has  been 
taken  (assuming  p  <<  1).  Comparison  of  (8)  and  (9)  shows 
that 


dP 


^ 


Use  of  this  approximation  to  the  derivative  gave  much 
improved  convergence  with  the  above  starting  value.  The 
routine  (POFKINV)  is  quite  simple  and  converges  well 
over  a  broad  range,  even  when  the  condition  for  (9)  does 
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not  hold.  In  less  than  15  iterations  it  will  match  the 
input  value  of  P(>k)  to  at  least  5  digits  relative 
accuracy.  No  trouble  was  experienced,  but  if  any  should 
be,  (9)  could  be  replaced  with  the  term  actually  biggest, 
which  is  just  (9)  with  k  replaced  by  the  smallest  integer 
^  p(M  +  1)-1,  which  might  be  done  recursively.  The  deriva¬ 
tive  might  also  need  to  be  changed  in  that  case. 

OTHER  CALCULATIONS:  Routines  were  written  to  compute  equa¬ 
tion  (2)  of  the  text  and  its  inverse,  but  they  are  not 
included  here  because  of  the  negative  results  and  their 
doubtful  utility.  The  inverse  routine  was  not  very  satis¬ 
factory  and  frequently  failed  to  converge,  although  more 
than  enough  highly  accurate  results  were  obtained  to 
reach  the  conclusions  stated  in  the  text. 

A  special  routine  was  written  to  compute  the  prob¬ 
ability  of  at  least  k  events  out  of  two  groups  having 
different  probabilities,  although  the  more  general  pro¬ 
cedure  of  Wirth  (1971)  could  also  have  been  used.  In 
the  notation  of  (2) 

f(;>o  *  l  pc”1|p1j 

or 

”l  M  .  M.-i 

P(ik)  •  l  q1)  P;  (1-p,)  1  S, 
i3i  L 
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where 


M 


S9  -< 


2  Mo  i 

l  C1Z)  P2  (i-P2) 
j=k-i  J  2  2 

1  ,  if  i>k 


M2-j 


,  if  i<k 


and  i  5  max(0,  k-M2).  The  routine  (P0FK3)  is  fairly 
straightforward.  False-alarm  rates  for  this  case  were 
computed  from  POFK  and  its  inverse,  since  false-alarm 
rates  were  assumed  to  be  the  same  for  all  stations,. 
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Function  qhcfc  f*el»m,n2  > 

ONCE  =  h'ROP<  )!>►  ).  WHt-PF  U  D1STHIH.  AS  NONl-UENTHAL  F  .  W/  Nor.- 
CEMTPALITY  PARA*.  ANO  *N1*N2*  UfcURfceS  OF  FRtEDOF.  (Nl  fVfcN), 

EL?  s  t-L  *.S  T  pNl  ■  Nl  |  s  N2 

em  *  e n «;  *.5  *  x  ■  en27<en2  ♦  eiu*F) 

Y»1.“X  %  jjH  I  I  I  |«*CH 

F I  *  1*  *  NUT  ■  Nl / 2  -  1 

Do  10  I  s  l.NIJT 

T  =  T  *  V  *  FM/Fi  F  SM  ■  SM  ♦  J  I  EM  «  fcM  ♦  1. 

10  F I  ■  F I  ♦  ]  . 

P  =  s*  *  AJ  «  fcL*  »  D  ■  1. 

DO  20  J  ■  1 >  ?00n 

T  r  T  *  Y  *  FM/F  I  1  *M  S  bM  ♦  1 

Tj  «  AJ  *  SM  %  p  s  p  ♦  jj 

IF (  T*/H.lt  .IE-7  )  Jo, 15 

15  Em  ■  EM  ♦  i .  »  Fl  ■  FI  ♦  1.  I  y  ■  d  ♦  i  ! 

20  AJ  «  AJ  *  EL?/T 
PRINT  1 

1  F  OPMA I ( >  2?H  pooh  CUN  VI-  PGf  HCE  UNCF) 

30  QNTF  ■  M*rvpF(-EL2> 

Return 

End 
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rUMCTlO'X  fjhCF  (  F  » EL » '  1>N?  »  Ml,? 

Data  <TP|  =  .  *-»6oiV  772M ) 

QMrf  =  PROP!  IJ ).  WFIFRF  u  pISTKlB.  AS  NON-UENTHAL  F,  W/  NQN- 
CewTRALlTY  PARA''.  *tl  *  ANn  *N1»N2*  UfcliRhcS  OF  FREEDOM,  |  N 1 . 2  Q  n  D  )  . 


El?  s  H-  *.s  i  gNl  ■  Nl 

X  =  FNP'IENP  ♦  CN1*F) 

T.'l  =  S  n  R  T  F  I  RN1*F/EN?  ) 

T  :  Sf  !  1,  %  pi  a  3. 


DO  10  I  *  1 ,  Mi  IT 
T  s  T  *  X  *  ( F  J-1.  >/F  I 
1 0  F  I  a  F  l  *  ? . 

A  s  TFI  *  (TWFT  *  X*TN*SM» 

G  z  TFI  *  T  *  ( F  l  -1  .  >  *  x#X*TN 

T  s  SF  3  i  .  P  F 1  *  3 . 

DO  20  i  *1 .MUT 

T  =  T  *  Y  *  fm/FI  %  qn  ■  bM  ♦  1 
20  Fl  a  F 1  ♦  ?. 

P  *  SM  ♦  AJ  «  tL2 


DO  30  J  «  1,2000 
T  a  T  *  Y  *  FM/F 1 
Tj  ■  AJ  *  SM 
IF (  TJ/P.LT.iF-7  ) 
25  Em  «  EM  ♦  ?. 


%  SM  *  SM  ♦  ) 

F  P  !  M  T,| 

*  F 1  *  F I  ♦  2 . 


30  Aj  *  AJ  *  FL?/n 
PR  I  A  T  1 


1  FORMaK/  ??M  PPCM  COFVFRGfNCE  CINCF) 

40  QNCF  «  A  ♦  F  X  P  F  <"EL2)*F*P 

Rftdra 

Ehp 


mw 

%  6N2  a  N2 

*  Y  «  1.-  X 

*  TPtfc  T  =  ATANF{1./TNJ 

1  NUT  a  (N2-3)/2 

*  SM  a  SM  ♦  I 


*  fcM  a  P  M2  ♦  1. 

*  NUT  a  <Nt-3>/2 

%  pM  a  p  m  ♦  2 . 

4  D  ■  1  • 


»  D  *  D  ♦  1  . 
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FUNCTION  OMCF(  F»EL»M.N?  >  N<i  hVr-N 

QMrF  s  K?np(  U>F  )#  u'Hf  RE  0  nlSTKIB.  AS  i\ON-CbNTNAL  F  >  W/  NON-  2 

CENTRALITY  Para''*  *fel  *  AND  *N1»Ni2*  y  fcORfccS  OF  FREEDOM,  (V  £  7  f  *  >  •  i 

Em  *  N1  *  FN?  *  N2  9 

X  s  F  N  2  /  ( F  M  ?  ♦  b  N 1  *  F  )  i  Y  «  1*-  X  o. 

EM4  *  CNI1  -  4.  1  FNI  V  ■  £N14  ♦  fcL*f  7, 

T  ?  *  0.  1  $m  *  T1  a  1.  %  F21  =  2.  n, 

nut  ■  N2/?  -i  v 

Do  10  I  a  l.MUT  ’  u 

T  S  »<  1 1  1.1 

T 1  ■  ((tNLY*?.*F2I >*Tl  -  ( fcNl44F  2  I  ) ***T  2  >  *  X/F21  2 
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FUMCT  1  ON  f  I  M  V  (  U»M»  •>  2  ) 

Q  =  FF0b(  '>FIfV  ),  i  mCTHld.  AS  F*  */  M.N2  DEGREES  OF  FREtDOM. 
lOVtFSt  (TUmilaM  OF  F-ni^THU..  tilVeb  T HnfcShOLD  AS  A  FUNC.  OF 
Fa!  Sf*aLi\Pi-  P4Tc  FOR  FtRwpFi  DETECTOR.  (ALL  N1.N2)  MW 


EM  e  Mi 

+ 

FN2  *  m2 

* 

R  s 

FM1/EN2 

Em  a  (EM  aF  N?  >*  .5 

T 

F- N 1 H  a  bNl*.5 

iFf  XI*  DUF  ( r.  1  ,  ?  )  1 

2  #  1 

1 

D  r  PN2  * .P 

• 

a— ♦ 

n 

'll 

it 

t 

NUT 

s  Nl/2 

QO  M  4 

2 

1 F  (  XMOuf  ;  f|?  ,  2  )  > 

6  • 

3 

D  *  FM  *  .5 

L  =  B  *  1J 

» 

NUT 

a  N2/2 

4 

DO  5  l  *  1 , NUT 

B  r  p  *  'I/O 

* 

ii  !  u  M> 

5 

D  =  I  ♦  i  . 

B  =  P  /  [) 

4 

(?D  TO  10 

6 

8*3  •i4*1  59?x,5.ffc 

V 

tl  *  .5 

% 

D  = 

1  • 

IT  *  ( Nx-1 ) /? 

► 

\UI  a  (IM14N2I/2  - 

1 

DO  S  I  *  i  ,NIT 
B  s  e  *  ii/p 
if (  i.eu.it  >  7» » 
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Function  quantf‘p>  id 

P  ■  PR  OB  (  IJ<0U*MF  »#  WHERk  U  DISsTRlB.  AS  ZEHQ-MfcAN  UNIT  NORMAL.  2() 

Inverse  of  prorf  mw  jo 

SigN*1.I0*P  4  (j 

IFIP.LE*  n.  .OR  P.UE.l.)  If), 20  bo 

10  QUANTF  «  P  60 

PRINT  1»  P  J,  RETURN  7  0 

1  FORMA  T ( 2<h  PRRrR  IN  nUANTF<PI#  Pa  fc20.10)  8p 

20  1F(  P.Gl  .  n  .5  j  25,30  Vu 

25  SIGN  ■  -t  .  %  0  «  1 P  IOC 

30  T  ■  SCRIM-?.*  LOaF(QJ)  Ilf 

OUANTF  ■  S!GN*(  <2.&l5*>i7  ♦  T*(.«Q?a?J  ♦  T*  .  U10328 )  )/ 1 1 ,  ♦  T*  1 2 1. 

*  (1 • 43?7flfl  ♦  T*<.lP9?6*  ♦  I*.U01Auem  -  T  )  lJf 

RETURN  i4o 

End  l^d 
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function  pofk<  n#a,p  j 

Sum  of  dlNOMUl  UJblwlH. 

P"FK  *  PROP.  OF  *T  LEAST  K  EVENTS  UUT  CF 
wHfN  ALL  FVENTS  HAVE  PRflfl.  P 
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POFA  *  I  ■  p  •  *  f. 
R  «  < l . -p  »  /  p 
KP  i  M  1 
50  IP  I  i  KP.N 
T  *  T  *  U/B)  * 
A  «  A  -  j  . 
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END 
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*0 
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Ho 
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End 
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FUNCTION  P()FK3<  K*Nl*N?,  Pl*P2  ) 
pBOb.  OF  »K  TOTAL  OF  Twn  StTS*  TH 
W/  PFCH*  Pi*  TNfc  SbCCND  KJ>  bL£*E* 
N*  l  M  AM)  K  S  iU«N? 


F  I HS  T  CONTAINING  Nj  FLFPFMP 
S  N/  M-00.  P2 


T 1  ■  FI  **  Nl  » 

R 1  ■  « 1 • -PI )/P1  ♦ 

plF«3  *  «*. 

A  ■  1 1  ■  N 1  ♦ 

JS  ■  t 

]F(  K.GI.N?  )  i«2 
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a  1  • 
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